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INTRODUCTION 

The Importance of aromatic hydrocarbons In today's practical hydrocarbon 
fuels is well known. This fact has resulted In several new studies of the 
oxidation and pyrolysis mechanisms of these compounds. A recent review paper 1 
has summarized the qualitative knowledge of the mechanism of oxidation of the 
simplest aromatics, benzene, toluene, and ethyl benzene. In addition, many of 
the individual reactions in the pyrolysis and oxidation of benzene and toluene 
have been studied in several papers. 2-8 This information has been used^ to 
construct the first detailed mechanism of benzene oxidation. The mechanism 
was tested using available experimental data on ignition delay times 10 and on 
temperature and composition profiles measured during benzene oxidation In a 
highly turbulent reactor. 1 It was found to compute the measured Ignition 
delay times and temperature profile reasonably wel 1 and to semiquantitatively 
reproduce some of the composition profiles. 

In the present work we develop a detailed mechanism for toluene oxidation 
using the same method as applied previously to the benzene oxidation. Recent 
Information^-^ on toluene pyrolysis and oxidation reactions was combined with 
the detailed benzene oxidation mechanism. The resulting mechanism was used to 
compute experimentally measured ignition-delay times 10 for shock-heated 
toluene-oxygen-argon mixtures and composition profiles for two toluene-oxygen- 
nitrogen mixtures in a turbulent flow reactor. 3 Most of the rate coefficients 
for the toluene reactions were used at their published literature values. 

Only those rate constants with large uncertainties were adjusted. The 


reactions controlling the ignition process and the profiles of various species 
concentrations were determined by an extensive sensitivity analysis using the 
new NASA Lewis Research Center chemical kinetics and sensitivity analysis 
code .n*12 

In the sections that follow we present comparisons of computed and 
experimental results and describe the sensitivity analysis results. 

TOLUENE OXIDATION MECHANISM 


The qualitative toluene oxidation and pyrolysis paths already outlined 1 
have been used along with the results of recent experimental work 3-6 to write 
the following initiation and chain propagation scheme Involving toluene and 


its pyrolysis fragments: 

Initiation: 

C7H8 + O2 

C?H 8 

C7 h 8 

Chain propagation: 

H + C 7 H 8 
H + C 7 H 8 

H + C7H 8 

ch 3 + C 7 H 8 
CH3 + C7H 8 

c 6 h 4 ch 3 

c 6 h 4 ch 3 

c 6 h 5 ch 2 


c 6 h 5 ch 2 + ho 2 

(1) 

C 6 H 5 CH 2 + H 

(2) 

c 6 h 5 + ch 3 

( 3 ) 

c 6 H 6 + ch 3 

( 4 ) 

c 6 h 5 ch 2 + h 2 
benzyl 

( 5 ) 

C 6 H 4 CH 3 + h 2 
methyl 

(6) 

phenyl 


c 6 h 5 ch 2 + ch 4 

(7) 

C6H 4 CH 3 + CH 4 

(8) 

c 4 h 3 + C 3 H 4 

( 9 ) 

C 3 H 3 + 2 C 2 H 2 

(10) 

C 3 H 3 + 2 C 2 H 2 

(11) 


2 


ch 3 + c 6 h 6 * ch 4 + c 6 h 5 


( 12 ) 


C3H4 + M C3H3 + H + M 

( 13 ) 

H + C3H4 C3H3 + H2 

( 14 ) 

0 + C 7 H 8 ^ OHC7H7 
cresols 

( 15 ) 

OH + C7Hg C6H5CH2 + H2O 

( 16 ) 

C 6 H 5 CH2 + 0 w* C6H5CHO + H 

( 17 ) 

C6H5CH2 + 02^ C7H7O + 0 
benzoxy 

( 18 ) 

OH + C7H7OH * C7H7O + H2O 

( 19 ) 

C 6 H 5 CH 2 + OH + M * C7H7OH + M 

benzyl 

alcohol 

(20) 

C5H5CH2 + HO2 ^ C7H7O + OH 

(21) 

C7H7O * C6H5CHO + H 

(22) 

C 6 H 5 CHO * C6H5 + HCO 

( 23 ) 

2 C 6 H 5 CH 2 ^ Ci 4 H 14 
bibenzyl 

( 24 ) 

OHC7H7 + H 2 ^ C6H5OH + ch 4 

( 25 ) 

2 C5H5 C12H10 

( 26 ) 

C 6 H 5 CHO + 0 2 ^ C6H5CO + H 0 2 

( 27 ) 

C 6 H 5 CHO + OH 5* C6H5CO + H2O 

( 28 ) 

C 6 H 5 CHO + H 0 2 C6H5CO + h 2 o 2 

( 29 ) 

c 6 h 5 co^ c 6 h 5 + CO 

( 30 ) 


The species OHC7H7 represents a composite mixture of ortho, para, and meta 
cresols. The reactions above were combined with the benzene oxidation and 
combustion mechanism developed previously^ to give a system of 143 reactions 
among 46 species. As In the case of benzene oxidation, one reaction from the 
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hydrogen-oxygen system was found to be Important In the toluene oxidation 
mechanism, namely: 

H + 0 2 ^ OH + 0 (47) 

This reaction Is Important In both Initiation and chain propagation. A 
listing of the significant toluene and benzene Initiation and chain reactions 
(plus the H + 02 ^ OH + 0 reaction) Is given In Table I along with the rate 
coefficients used. All other reactions used In the computations have been 
given In our previous work. 9 Rate coefficients for the pyrolysis reactions 3, 
4, and 6 to 14 were all taken unchanged from the work of Pamldlmukkala et. al. 4 
Selection of the rate coefficients used for reactions 2, 5, and the oxidation 
reactions Is described below. 

SENSITIVITY ANALYSIS STUDY 

We used the method described previously 9 to compute normalized 
sensitivity coefficients of species concentrations, temperature, and pressure 
with respect to the parameters of the modified Arrhenius rate coefficient 
expression, k * AT n exp(-E/RT), for each reaction. Sensitivities with respect 
to the preexponentlal factors, A, of several benzene and toluene reactions are 
listed In Tables II and III. It should be noted that sensitivity coefficients 
with respect to the activation energy, E, are approximately equal In sign and 
magnitude to those for A under all conditions studied. Ignition delay 
times 10 were measured from pressure-time traces by a method described below. 
Shown In Table II are pressure sensitivity coefficients for Ignitions of two 
shock-heated toluene-oxygen-argon mixtures. These are the lowest temperature 
lean mixture and a high temperature stoichiometric mixture. As shown In the 
table the reaction of toluene with molecular oxygen (reaction 1) and that of 
hydroperoxyl radical with benzyl radical (reaction 21) are the most Important 
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reactions which control the ignition delay time. Also showing significant 
sensitivity here are the H + O 2 chain branching reaction (reaction 47), oxygen 
atom plus benzyl radical (reaction 17), and two Important benzene oxidation 
reactions of phenyl and cyclopentadienyl radical by molecular oxygen 
(reactions 38 and 45). The latter are of equal or greater Importance than the 
oxygen atom plus benzyl reaction in both of the mixtures shown in Table II. 

The direct reaction of benzyl radical with molecular oxygen Is insignificant, 
in direct contrast to the dominance of the phenyl plus oxygen reaction In the 
oxidation of benzene 9 caused by its much higher rate coefficient and molar 
rate. The conclusions drawn from Table II agree closely with the findings of 
previous investigators 3 that radical-radical reactions should be Important in 
the toluene oxidation and that the low rate of the benzyl-molecular oxygen 
reaction makes it unimportant. We note that reaction 5, though unimportant for 
the low temperature mixture, becomes as sensitive as reaction 17 in the higher 
temperature stoichiometric mixture. 

Table III presents sensitivity coefficients of several species concentra- 
tions for the oxidation of a lean toluene mixture in a turbulent reactor. 

These results show that the five reactions Important for ignition delay times 
are also Important in determining the concentration profiles of several 
species. In addition, reactions 4, 5, 15, 16, 20, and 29 have strong effects 
on the cresols, benzene, phenol, benzyl alcohol, and bibenzyl profiles. The 
phenyl and cyclopentadienyl oxidation reactions 38 and 45 have a moderate 
effect on the concentration profiles of phenol and carbon monoxide. As only 
an estimate of the rate coefficient of the highly sensitive reaction between 
toluene and molecular oxygen exists, 8 we consider this reaction to be an 
adjustable parameter to be used in matching the computations to the 
experimental results. 
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The following procedure was used to obtain the best possible matching to 
the experimental data: Reaction 21 was set within 20 percent of its collision 

theory estimate. 3 The activation energy and preexponential factor for reac- 
tion 1 were adjusted to closely predict the high-temperature (1343 to 1600 K) 
experimental ignition delay times in argon-diluted mixtures and also give a 
temperature rise of no more than 5 K for the two different mixtures in the 
turbulent flow reactor at 1180 K. In these computations the rate coefficient 
of the pyrolysis reaction 2 was taken as one half that of reaction 3. This is 
the ratio given by Pamidimukkala et al. 4 for their lowest temperature of 
1600 K and was used without any attempted variation. The rate constant of 
Robaugh and TsangS was used for reaction 5 because it is a directly measured 
experimental value and gave the best compromise in attempting to match the 
highest temperature Ignition delay times and some of the turbulent reactor 
experimental composition profiles. Finally, the rate coefficients of reactions 
15 to 30 were either estimated or used at or near their literature estimates to 
give the best possible agreement with the turbulent flow reactor data. Only 
those coefficients which had no effect on the computed pressure profiles were 
changed. 


DESCRIPTION OF COMPUTATIONAL PROCEDURE 

The same procedures used in our benzene mechanism study^ were used in 
this work to model the experimental ignition delay time measurements 10 and the 
experimental composition profile data 3 that have been published for toluene- 
oxygen mixtures. Only the important details will be summarized here. 
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Shock Tube Ignition Experiments 


A constant volume batch reaction model was used for the shock-heated 
mixtures. The reported Initial reflected-shock temperature and pressure 
conditions 10 were recomputed, as described previously, 9 applying a small 
correction for attenuation of the shock velocity to each data point. As shown 
by Brabbs and Robertson, 13 all data points with Ignition delay times less than 
100 psec were considered Inaccurate and eliminated from consideration. 
Experimental Ignition delay time was determined from each experimental 
pressure versus time curve as the time of the first "significant" rise In the 
pressure. 10 Each computed time was measured from the corresponding pressure 
versus time plot as described previously. 9 The Ignition delay time 
represented a pressure rise of about 3 to 6 percent over the Initial value. 

The thermodynamic data used for all computations are from the NASA Lewis data 
base, which Is part of the Gordon and McBride Chemical Equilibrium Code. 14 
New, Improved data for many aromatic species were kindly provided by Bonnie J. 
McBride of this laboratory. 

Turbulent Flow Reactor 

As described previously, 9 the turbulent flow reactor was modeled as a 
constant pressure homogeneous batch reaction. A detailed description of the 
reactor Is given by Hautman, 15 who Indicates that the reactor was run at a 
constant pressure of 1 atm. In this apparatus fuel Is Injected Into a 
nitrogen-diluted, highly turbulent stream of oxygen. The exact zero of 
reaction time Is unknown and was taken as the point of fuel Injection Into the 
hot oxidant stream. Distance profiles were converted to time profiles by use 
of the measured flow velocities In the reactor. 
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RESULTS AND DISCUSSION 


Comparison of Computed and Experimental Ignition Delay Times 

Comparisons of experimental and computed ignition delay times are shown 
for four different starting mixtures in Figs. 1 to 4, where logarithm of 
ignition delay time is plotted against the reciprocal of temperature. Initial 
conditions are given in Table IV along with a comparison of computed and 
experimental results for all data points. Included in Table IV is an error 
analysis of the results. The percent difference between each experimental and 
computed ignition delay time is given along with the percent standard deviation 
defined previously 9 for each of the four sets of experimental conditions. 

Figures 1 to 4 show fair to good agreement between computed and 
experimental results, better for the stoichiometric mixtures than for the lean 
mixture. Shown are the individual computed and experimental points as well as 
least-squares lines for each set of points, fitted to the empirical equation 

x » A e^/RT 

where x is the ignition delay time (experimental or computed), R is the 
universal gas constant and AE is the activation energy term for each set of 
initial conditions. For mixtures 1, 2, and 4, with initial pressures around 
2 atm, the computed temperature dependence is weaker than that observed 
experimentally. For mixture 3, the dilute stoichiometric mixture with initial 
pressures around 6 atm, the computed temperature dependence is stronger than 
that observed experimentally. Computed and experimental activation energies 
are tabulated in Table V. The experimental activation energy for mixture 3 is 
significantly lower than that for mixture 2, whereas the computed activation 
energy for mixture 3 is only slightly lower than the value for mixture 2. 
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This experimental result Is quite different from that for benzene-oxygen 
Ignitions, 9 which had very similar Initial shock conditions to those for the 
toluene experiments. For benzene neither the experimental nor the computed 
temperature dependences changed with Initial pressure. In fact both the 
experimental and computed benzene-oxygen activation energies changed only 
moderately for all four Initial mixtures used. 

A comparison of Ignition delay time measurements for mixtures 2 and 4- 
shows the effect of argon dilution for a constant equivalence ratio of 1.0. 
Figure 5 shows computed and experimental results for these two mixtures. Only 
the least-squares lines from Figs. 2 and 4 are shown for clarity. Our computed 
results can be seen to satisfactorily match the magnitude of the experimentally 
observed effect of argon dilution. In our previous benzene mixture 
computations ,9 the computed effect of argon dilution was much smaller than the 
experimentally observed effect. 

A comparison of results for mixtures 2 and 3 shows the effect of 
Increasing initial molar concentrations by shock-heating of the same molar 
mixture at two different initial pressures, 2 and 6 atm. The least-squares 
lines of Figs. 2 and 3 are replotted in Fig. 6, which shows that the computed 
magnitude of this concentration effect matches the experimental magnitude 
better at low temperatures than at high temperatures. 

In summary, these comparisons have shown that our proposed toluene 
oxidation mechanism reasonably matches the experimental ignition delay time 
data over a wide range of initial conditions. The agreement between 
computation and experiment is about the same as that obtained for benzene 
mixtures in our previous study. 9 
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Comparison of Computed and Experimental Turbulent Reactor Results 

The turbulent-reactor toluene oxidation experiments were performed at 
essentially constant temperature, and no temperature versus time profile was 
reported. 3 The maximum measured temperature rise was reported to be 5 K in 
all experiments. As stated previously, we adjusted the rate coefficient 
parameters of reaction 1 so that the computed maximum temperature rise at the 
outlet of the reactor agreed with the experimental value. 

Figures 7 to 14 show computed and experimental composition versus time 
profiles for toluene oxidation in the flow reactor with initial temperature of 
1180 K and a pressure of 1 atm. Two mixtures, with equivalence ratios (cp) of 
<p * 0.63 and 1.4, were used in the experimental study and each mixture 
contained 0.14 mol % toluene. Figure 7 shows toluene versus time profiles for 
both mixtures. The computed profile for cp = 0.63 matches the experimental 
profile well, except toward the end of the reaction. Reasonable agreement is 
also obtained for the <p = 1.4 mixture, for which the maximum difference is 
about 15 percent. For both mixtures there is satisfactory agreement between 
experimental and computed slopes of the curves at early reaction times. This 
prediction of the fuel versus time profile contrasts with the failure to 
accurately predict a benzene destruction profile measured In the same reactor. 
Our benzene oxidation mechanism 9 predicts much more rapid destruction of 
benzene than is observed experimentally. Prediction of other concentration 
profiles was not as successful. The computed phenol profile (Fig. 9) for 
cp * 1.4 gives fair quantitative matching to the experimental profile. However, 
the computed profile for cp = 0.63 and the other computed profiles for cresols 
(Fig. 8), benzene (Fig. 10), benzaldehyde (Fig. 11), carbon monoxide (Fig. 12), 
benzyl alcohol (Fig. 13), and bibenzyl (Fig. 14) show only qual itative 
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agreement with the experimental profiles. This qualitative agreement has been 
obtained even though there is little experimental information available about 
the reactions of many of these species. Several attempts were made to improve 
the overall agreement of the species profiles by certain rate coefficient 
variations. However, it was found that the single changes improved some 
profiles and made others worse. The results given appear to be the best 
compromise using the given set of reactions. 

In summary, our proposed mechanism reasonably predicts the destruction 
profile of toluene, but is only partially successful in matching other 
experimental concentration profiles. 

CONCLUDING REMARKS 

We have presented a detailed toluene oxidation mechanism which reasonably 
computed measured ignition delay times in argon-diluted mixtures over a wide 
range of experimental conditions. In addition, the mechanism computed fairly 
good toluene versus time concentration profiles for the nitrogen-diluted 
oxidation in a turbulent flow reactor. Profiles of several other species 
concentrations were qualitatively matched. 

Sensitivity analysis shows that the direct reaction of toluene with 
molecular oxygen strongly effects the profiles of temperature, pressure, and 
many species concentrations. This is in sharp contrast to the unimportance of 
the corresponding reaction in the benzene oxidation. A comparison of the heat 
of reaction for the benzene plus oxygen reaction (~60 kcal/mol) with that for 
the toluene plus oxygen reaction <~35 kcal/mol) justifies the much lower 
activation energy and higher rate coefficients for the latter reaction. 
Computations show that the molar rate of the toluene-oxygen reaction is always 
several orders of magnitude greater than that of the benzene-oxygen reaction 
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for similar temperature and molar concentration conditions. The much higher 
reactivity of toluene with oxygen accounts for this reaction's being a major 
path for toluene oxidation. The benzyl-molecular oxygen reaction was found to 
be quite unimportant in toluene oxidation. This contrasts with the dominant 
effect of the corresponding phenyl-oxygen reaction in the oxidation of benzene. 
These facts are consistent with the idea that benzyl radical is conjugatively 
stabilized 1 and is less reactive than phenyl radical. The benzyl reaction with 
molecular oxygen is endothermic, whereas the corresponding phenyl reaction is 
very exothermic. The latter reaction, with its much higher rate coefficient 
and molecular rate, is one of the most important reactions in benzene 
oxidation, whereas the benzyl-oxygen reaction has very little effect on toluene 
oxidation. Because benzyl is a relatively stable radical, its reactions with 
other radicals, primarily hydroperoxyl , are its Important ones in the oxidation 
of toluene. 

A rate coefficient expression for the toluene-molecular oxygen reaction 
was found which predicted the temperature dependence of ignition delay times at 
high temperature (1300 to 1600 K) and also matched the very small temperature 
rise reported for the turbulent reactor at 1180 K. The results of this study 
have given a toluene oxidation mechanism that can be used for ignition and 
combustion modeling in practical, well-mixed combustion systems. 
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TABLE I. - TOULENE OXIDATION MECHANISM 


Num- 

ber 

Reaction 

A 

CGS units 3 

n 

— 

E 

cal / mol 

Reference 

1 

C 7 H 8 + 0 2 ^ c 6 h 5 ch 2 + ho 2 

3 -30x1 o ' 4 

0.0 

38 000 

Thi s 

work 

2 

c 7 h 8 ^ C 6 H 5 CH 2 + H 

4.45X10 12 



72 600 

4 

3 

c 7 h 8 ^ c 6 h 5 + ch 3 

8. 91 xlO 12 



72 600 

4 

4 

H + C 7 H 8 ^ C 6 H 6 ♦ CH 3 

4-OOxlO 13 



5 120 

4 

5 

H + C 7 H 8 ^ C 6 H 5 CH 2 + H 2 

1 . 20 xl 0' 4 



8 220 

5 

6 

H + C 7 H 8 ^ C 6 H 4 CH 3 + H 2 

2.51 xlO 14 



16 000 

4 

7 

CH 3 + C 7 H 8 ^ C 5 H 5 CH 2 + CH 4 

4. 37x1 O - 4 

5.0 

8 300 



8 

ch 3 + c 7 h 8 ^ c 6 h 4 ch 3 f ch 4 

4.37X10 ' 4 

5.0 

12 300 



9 

c 6 h 4 ch 3 ^ c 4 h 3 + c 3 h 4 

1 .OOxlO 16 

0.0 

82 000 



10 

C 6 H 4 CH 3 ^ C 3 H 3 + 2 C 2 H 2 

1 .OOxlO 16 

0.0 

83 000 



11 

C 6 H 5 CH 2 ^ C 3 H 3 + 2 C 2 H 2 

1 -78xl0 14 

0.0 

84 800 



12 

CH 3 + C 5 H 5 ^ CH 4 + CgH 5 

4.37xlO - 4 


5.0 

12 300 



13 

M + C 3 H 4 ^ C 3 H 3 + H 

2. OOxlO 17 


).0 

65 000 



14 

H + C 3 H 4 ^ C 3 H 3 + H 2 

6.92xl0 14 



14 500 

J 


15 

0 + C 7 H 8 ^ OHC 7 H 7 

2 . 20 xl 0 13 



3 800 

Adj from 7 

16 

OH + C 7 H 8 ^ C 5 H 5 CH 2 + H 2 0 

3. OOxlO 12 



3 000 

Adj from 8 

17 

C 6 H 5 CH 2 + 0 -C 6 H 5 CHO + H 

T . 20 x 10' 3 




Adj from 3 

18 

C 6 H 5 CH 2 + O 2 ^€7870 + 0 

6.30X10 12 



43 000 

3 

19 

OH + C 7 H 7 OH ^ C 7 H 7 O + H 2 0 

1 -OOxlO 13 



5 000 

Estimated 

20 

C 6 H 5 CH 2 + OH + M -C 7 H 7 OH + M 

1. OOxlO 17 




Estimated 

21 

C 6 H 5 CH 2 + H0 2 ^ C 7 H 7 O + OH 

3.60X10 12 




Adj from 3 

22 

C 7 H 7 O -CgHsCHO + H 

1 .OOxlO 12 






23 

C 6 H 5 CH0 ^ C 6 H 5 + HCO 

1 .OOxlO 16 



82 000 

Estimated 


a Parameter$ in the expression k = AT n exp(-E/RT). 
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■JFjS 

w 

Reaction 

A 

CGS units 

n 

E 

cal /mol 

Reference 

24 

2 C 6 H 5 CH 2 ^ C 14 H 14 

1 .OOxlO 14 

0.0 


Estimated 

25 

OHC 7 H 7 + h 2 ^ c 6 h 5 oh + CH 4 

4. OOxlO 11 



5 000 



26 

2 C 5 H 5 ^ C 12 H 10 

1 .OOxlO 12 






27 

c 6 h 5 cho + o 2 ** c 6 h 5 co + ho 2 

5. OOxlO 12 



35 000 



28 

c 6 h 5 cho + OH * c 6 h 5 co + h 2 o 

5. OOxlO 12 



5 000 



29 

c 6 h 5 cho + ho 2 ** c 6 h 5 co + h 2 o 2 

1 .OOxlO 14 



5 000 



30 

c 6 h 5 co * c 6 h 5 + CO 

4. OOxlO 14 



29 400 

16 

31 

c 6 h 6 + 0 2 * C 6 H 5 ♦ H0 2 

6.31xl0 13 



60 000 

8 

32 

C 6 H 6 - C 6 H 5 + H 

5. OOxlO 15 



108 000 

17 

33 

C 6 H 6 + H ^C 6 H 5 + H 2 

2 . 50x 1 0 1 4 



16 000 

2 

34 | 

C 6 H 6 + 0 - C 6 H 5 + OH 

2 . 78X10 1 3 



4 910 

7 

35 

CgHg + OH ^ C 6 H 5 + H 2 0 

2.13xl0 13 



4 580 

18 

36 

C 4 H 3 + M ^ C 4 H 2 + H + M 

3.31xl0 51 

- 10.0 

63 000 

2 

37 

CgHsO ^ C 5 H 5 + CO 

2 .51 xlO 1 1 

0.0 

43 900 

19 

38 

c 6 H 5 + o 2 ^ c 6 h 5 o + 0 

4.50X10 12 



15 000 

9 

39 

C 6 H 5 ^ C 4 H 3 + c 2 h 2 

1 .58xl0 15 



82 000 

2 

40 

C 6 H 5 0H ^ C 5 H 5 O + H 

6 . OOxlO 13 



88 000 

9 

41 

! C 5 H 6 + 0 ^ C 5 H 5 0 + H 

5. OOxlO 12 



10 000 



42 

C 6 H 5 OH + OH ** C 6 H 5 0 + H 2 0 

8 . OOxlO 12 



5 000 



43 

CgHs + CgHg ^ C6H5 + C5H6 

2. OOxlO 11 


■ 

10 000 



44 

C5H5O + M ^ C4H5 + CO + M 

7.59X10 13 



9 



45 

C 5 H 5 + 0 2 ^ C 5 H 5 O + 0 

2. OOxlO 12 


■ 

' u 



46 

C 4 H 5 ^ c 2 h 3 + c 2 h 2 

1 ,40xl0 13 



32 .900 



47 

H + 0 2 ^ OH + 0 

1 . 66 xl 0 14 



16 400 

20 


15 











TABLE II. - PRESSURE SENSITIVITY COEFFICIENTS FOR SHOCK IGNITION OF 
TOLUENE-OXYGEN-ARGON MIXTURES 


Num- 

ber 

Reaction 

Normalized pressure sensitivity coefficient. 
A/p (3p/3A) 

<p = 0.331, T = 1334 K 

<p - 1.0, T - 1535 K 

1 

C 7 H 8 + 02 ^ C 6 H 5 CH 2 + H 02 

0.02506 

0.02423 

2 

C 7 H 8 - C 6 H 5 CH 2 + H 

-0.00258 

-0.00056 

3 

C 7 H 8 ^ C 6 H 5 + CH 3 

-0.00282 

0.00068 

4 

H + C 7 H 8 * C 6 H 6 + CH 3 

-0.00461 

-0.00391 

5 

H + C 7 H 8 ** C 6 H 5 CH 2 + H 2 

-0.00396 

-0.00763 

15 

0 + C 7 H 8 ^ OHC 7 H 7 

-0.00959 

-0.00659 

16 

OH + C 7 H 8 ^ C 6 H 5 CH 2 + H 2 0 

-0.00762 

-0.00493 

17 

C 6 H 5 CH 2 + 0 ^ C 5 H 5 CHO + H 

0.00842 

0.00732 

18 

C 6 H 5 CH 2 + 0 2 ^ C 7 H 7 O + 0 

0.00005 

0.00009 

20 

C 6 H 5 CH 2 + OH + M - C 7 H 7 OH + M 

-0.00107 

0.00013 

21 

C 6 H 5 CH 2 + H0 2 C 7 H 7 O + OH 

0.01733 

0.02058 

23 

C 5 H 5 CHO ^ C 5 H 5 + HCO 

0.00209 

0.00402 

29 

C 6 H 5 CHO + H0 2 ^ CgHsCO + H 2 0 2 

0.00601 

0.00106 

38 

C 6 H 5 + 0 2 ^ C 6 H 5 O + 0 

0.00734 

0.00974 

45 

C 5 H 5 + 0 2 ^ C 5 H 5 O + 0 

0.01251 

0.00869 

47 

H + 0 2 ^ OH + 0 

0.01057 

0.01286 




TABLE 1 1 1 ( a ) . - SENSITIVITY COEFFICIENTS OF FOUR SPECIES IN A TURBULENT FLOW 
REACTOR AT EQUIVALENCE RATIO 0.63 
[Temperature = 1180 K; reaction time = 55 msec. 3 


Num- 

ber 

Reaction 

Normalized sensitivity coefficient of 
species concentration c 
A/c (3c/8A) 

Toluene 

Cresol s 

Phenol 

Benzene 

1 

C 7 H 8 + 0 2 ^ C 6 H 5 CH 2 + HO 2 

-0.2143 

0.4324 

0.4815 

0.1362 

2 

C 7 H 8 - C 6 H 5 CH 2 + H 

0.0416 

0.2706 

-0.2041 

-0.3026 

3 

C 7 H 8 ^ C 5 H 5 + ch 3 

0.0129 

-0.0507 

-0.0574 

-0.0255 

4 

H ♦ C 7 H 8 ^ C 6 H 6 ♦ CH 3 

0.0342 

-0.2437 

-0.4312 

0.7074 

5 

H + C 7 H 8 ^ C 6 H 5 CH 2 + H 2 

-0.0082 

-0.4349 

0.5212 

-0.1585 

15 

0 + C 7 H 8 - ohc 7 h 7 

- 0.0101 

0.3538 

-0.1150 

-0.1512 

16 

OH + C 7 H 8 ^ C 6 H 5 CH 2 + H 2 0 

-0.0366 

-0.2223 

-0.0858 

0.1251 

17 

c 6 h 5 ch 2 + 0^ C 6 H 5 CH0 + H 

-0.0016 

-0.2145 

0.1175 

0,1552 

18 

C 6 H 5 CH 2 + 0 2 * C 7 H 7 O + 0 

- 0.0002 

0.0011 

0.0008 

0.0005 

20 

C 6 H 5 CH 2 + OH + M ^ C 7 H 7 OH + M 

0.0016 

0.0183 

0.0354 

0.0379 

21 

C 6 H 5 CH 2 + H0 2 ** C 7 H 7 O + OH 

-0.1438 

0.1514 

0.3644 

0.2152 

23 

C 6 H 5 CHO ^ C 6 H 5 + HCO 

-0.0016 

0.0543 

0.0906 

0.0755 

25 

OHC 7 H 7 + H 2 ** C 6 H 5 OH + CH 4 

0.0091 

-0.1617 

0.1740 

-0.0231 

29 

c 6 h 5 cho + H0 2 ^ c 6 h 5 co + h 2 o 2 

-0.0681 

0.2788 

0.0959 

-0.0876 

35 

c 6 h 6 + OH ^ c 6 h 5 + h 2 o 

-0.0131 

0.1782 

0.1280 

-0.3712 

38 

c 6 h 5 + o 2 ^ c 6 h 5 o + 0 

-0.0238 

0.1207 

0.1272 

0.0069 

45 

C 5 H 5 + 0 2 ^ C 5 H 5 O + 0 

-0.0253 

0.1377 

0.2049 

0.0719 

47 

H + 0 2 - OH + 0 

-0.0562 

0.3721 

-0.0152 

-0.1576 



TABLE 111(b) . - SENSITIVITY COEFFICIENTS OF FOUR SPECIES IN A TURBULENT FLOW 
REACTOR AT EQUIVALENCE RATIO 0.63 
[Temperature = 1180 K; reaction time = 55 msec.] 


Num- 

ber 

Reaction 

Normalized sensitivity coefficient of 
species concentration c 
A/c ( 3c/3A) 

Benzal- 

dehyde 

Carbon 

monoxide 

Benzyl 
al cohol 

Bi benzyl 

1 

C 7 H 8 + 0 2 ^ C 5 H 5 CH 2 + ho 2 

0.4180 

0.6536 

0.6280 

0.1349 

2 

C 7 H 8 ^ C 6 H 5 CH 2 + H 

-0.0890 

0.1193 

-0.0686 

-0.3255 

3 

C 7 H 8 ^ C 5 H 5 + CH 3 

-0.0141 

-0.0564 

-0.0074 

0.0252 

4 

H ♦ C 7 H 8 ^C 6 H 6 + CH 3 

-0.0752 

-0.2625 

-0.1130 

0.1202 

5 

H + C 7 H 8 ^ C 5 H 5 CH 2 + H 2 

0.0230 

-0.0271 

- 0.0010 

0.1926 

15 

0 + C 7 H 8 OHC 7 H 7 

-0.2460 

-0.1817 

-0.0761 

0.1424 

16 

OH + C 7 Hg ** C 6 H 5 CH 2 + H 2 0 

0.0769 

-0.1558 

-0.4884 

0.4777 

17 

CgH 5 CH 2 + 0 ^ C 6 H 5 CHO + H 

0.2867 

0.1470 

0.0461 

-0.1727 

18 

C 6 H 5 CH 2 + o 2 ^ c 7 h 7 o + 0 

0.0009 

0.0006 

0 . 0002 

-0.0007 

20 

C 6 H 5 CH 2 + OH + M ^ C 7 H 7 OH + M 

- 0.0212 

0.0136 

0.5145 

-0.1224 

21 

C 6 H 5 CH 2 + H0 2 ^ C 7 H 7 0 + OH 

0.8072 

0.3246 

0.3505 

-0.0163 

23 

CgHsCHO ^ C 5 H 5 + HCO 

-0.0792 

0.0778 

- 0.0102 

-0.0744 

25 

OHC 7 H 7 + H 2 ^ C 5 H 5 OH + ch 4 

-0.0144 

-0.0234 

-0.0138 

0.0212 

29 

CeH 5 CHO + H0 2 ^ CgH 5 CO + H 2 0 2 

-0.3922 

0.3123 

0.2812 

0.1544 

35 

C 6 H 6 + OH ^ C 6 H 5 + H 2 0 

0.0397 

0.1488 

-0.0413 

-0.1053 

38 

c 6 h 5 + 0 2 ^ C 6 H 5 0 + 0 

0.0324 

0.1268 

0.0215 

-0.0500 

45 

C 5 H 5 + 0 2 ^ C 5 H 5 O + 0 

0.0457 

0.2177 

0.0658 

-0.0767 

47 

H + 0 2 ^ OH + 0 

0.1125 

0.0804 

0.1464 

0.0278 | 
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IGNITION DELAY TIME, pSEC 


TABLE V. - COMPARISON OF COMPUTED AND EXPERIMENTAL ACTIVATION 
ENERGIES FOR TOLUENE-OXYGEN-ARGON IGNITION DELAY TIMES 


Mixture 

description 

Activation energy, cal /mol 

Percent 

difference 

Experimental 

Computed 

No . 1 : c p = 0.331 

P = 2 atm 

61850 

43010 

-30.4 

No. 2: <p = 1 .0, 95% Ar 

P s 2 atm 

61770 

53210 

-13.9 

No. 3: <p = 1 .0, 95% Ar 

P s 6 atm 

38790 

46940 

21.0 

No. 4: <p = 1 .0, 85% Ar 

P = 2 atm 

53260 

41020 

-23.0 



1930 1390 1350 1315 

TEMPERATURE, K 

FIGURE 1. - IGNITION DELAY TIME VER- 
SUS RECIPROCAL OF TEMPERATURE FOR 
TOLUENE-OXYGEN-ARGON; MIXTURE 1, 
EQUIVALENCE RATIO = 0.331, INITIAL 
PRESSURE 3 2 atm. 



1610 1515 1930 

TEMPERATURE, K 

FIGURE 2. - IGNITION DELAY TIME VERSUS 
RECIPROCAL OF TEMPERATURE FOR TOLUENE- 
OXYGEN-ARGON; MIXTURE 2: EQUIVALENCE 

RATIO = 1.0, 95 PERCENT Ar. INITIAL 
PRESSURE = 2 ATM. 
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IGNITION DELAY TINE, jjSEC IGNITION DELAY TIME, (1SEC 



1610 1515 1430 1350 

TEMPERATURE, K 

FIGURE 3. - IGNITION DELAY TIME VERSUS 
RECIPROCAL OF TEMPERATURE FOR TOLUENE- 
OXYGEN-ARGON; MIXTURE 3: EQUIVALENCE 

RATIO = 1, 95 PERCENT Ar, INITIAL 
PRESSURE = 6 ATM. 



1610 1515 1430 1350 

TEMPERATURE, K 

FIGURE 5. - IGNITION DELAY TIME VESUS 
RECIPROCAL OF TEMPERATURE FOR TOLUENE- 
OXYGEN-ARGON; EFFECT OF ARGON DILUTION 
FOR EQUIVALENCE RATIO - 1.0; MIXTURE 
2, 95 PERCENT Ar; MIXTURE 4, 85 PER- 
CENT Ar; INITIAL PRESSURE = 2 ATM. 



1515 1430 1350 

TEMPERATURE, K 

FIGURE 4. - IGNITION DELAY TIME VERSUS 
RECIPROCAL OF TEMPERATURE FOR TOLUENE- 
OXYGEN-ARGON; MIXTURE 4: EQUIVALENCE 
RATIO = 1.0, 85 PERCENT Ar INITIAL 
PRESSURE = 2 atm. 



1610 1515 1430 1350 

TEMPERATURE, K 

FIGURE 6. - IGNITION DELAY TIME VERSUS 
RECIPROCAL OF TEMPERATURE FOR TOLUENE- 
OXYGEN-ARGON; EFFECT OF INITIAL REAC- 
TANT MOLAR CONCENTRATION: EQUIVALENCE 

RATIO = 1.0. MIXTURE 2: 95 PERCENT 

Ar; INITIAL PRESSURE = 2 ATM. MIXTURE 
3: 95 PERCENT Ar; INITIAL PRESSURE 

= 6 ATM. 
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TOLUENE MOLE FRACTION x 10 



FIGURE 7. - TOLUENE VERSUS TIME PROFILES 
FOR TOLUENE-OXYGEN REACTION IN NITRO- 
GEN. p = 1 atm, T = 1180 K. 



TIME, MSEC 


FIGURE 9. - PHENOL VERSUS TIME PROFILES 
FOR TOLUENE-OXYGEN REACTION IN NITRO- 
GEN. p = 1 ATM, T * 1180 K. 



FIGURE 8. - CRESOLS VERSUS TIME PROFILES 
FOR TOLUENE-OXYGEN REACTION IN NITRO- 
GEN. P = 1 atm, T * 1180 K. 



TIME, MSEC 


FIGURE 10. - BENZENE VERSUS TIME PROFILES 
FOR TOLUENE-OXYGEN REACTION IN NITRO- 
GEN. p = 1 atm, T = 1180 K. 
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BENZALDEHYDE MOLE 

BENZYL ALCOHOL MOLE FRACTION x 10 5 FRACTION x 10 5 



TIME, MSEC 

FIGURE 11. - BENZALDEHYDE VERSUS TIME 
PROFILES FOR TOLUENE-OXYGEN REACTION 
IN NITROGEN, p = 1 ATM, T = 1180 K. 



0 AO 80 120 160 

TIME, msec 

FIGURE 13. - BENZYL ALCOHOL VERSUS TIME 
PROFILES FOR TOLUENE-OXYGEN REACTION 
IN NITROGEN, p = 1 ATM, T = 1180 K. 



0 AO 80 120 160 

TIME, MSEC 

FIGURE 12. - CARBON MONOXIDE VERSUS TIME 
PROFILES FOR TOLUENE-OXYGEN REACTION 
IN NITROGEN. P = 1 ATM, T = 1180 K. 



0 AO 80 120 160 


TIME, msec 

FIGURE 1A. - BIBENZYL VERSUS TIME PRO- 
FILES FOR TOLUENE-OXYGEN REACTION IN 
NITROGEN, p = 1 atm, T = 1180 K. 
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